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C^- ■ We consider the supergravity dual of the N = 1* theory at finite tempera- 

ture by applying the Polchinski-Strassler construction to the black D3 brane 



o 
o . 

solution of Type IIB supergravity. At finite temperature the 5-brane probe 
action is minimized when the probe falls to the horizon, although metastable 
minima with r » rn persist for a range of temperatures. Thermal effects 
on the 3-form source for the hypermultiplet mass m and its order m 2 back 
^S^ . reaction on the other fields of the IIB theory are computed. We find unique 

solutions which are regular at the horizon and have the correct behavior on 
the boundary. For fixed temperature T, the horizon shrinks for increasing 
$—i ' m 2 suggesting that there is a critical temperature separating the system into 

high and low temperature phases. In the high temperature phase 5-branes 
are unnecessary since there are no naked singularities. Using the order m 2 
correction to the horizon area we calculate the correction to the entropy to 
be AS = —0A71AN 2 m 2 T, which is less than the free field result. 
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I. INTRODUCTION 

The Maldacena conjecture relating classical supergravity to strongly coupled confor- 
mal gauge theories has had many successes 0. Perhaps more remarkable is the fact that 
one can apply it to nonconformal theories as well. The first, and perhaps most famous appli- 
cation has been to M = 4 super Yang-Mills at finite temperature |J . One may consider this 
as a gauge theory compactified on a circle where the fermions have anti-periodic boundary 
conditions, breaking the supersymmetry. This dimensionally reduced gauge theory is ex- 
pected to be confining with a mass gap. By arguing that this theory has a supergravity dual 
which is a black hole in a curved space, Witten was able to demonstrate both the existence 
of an area law and a mass gap in the strong coupling limit. 

One would like to find other models that exhibit confining behavior. Such examples have 
been found, either by adding spin to the black hole, breaking some of the global symmetries 
of the gauge theories [§]-§|, or by considering more exotic theories such as Type string 
theories [|9| [3~2|| or even type II theories with a dilaton turned on ||13HT5| , breaking the M = 4 



supersymmetry. In these latter cases, while there appear to be gaps and area laws, there are 
also naked singularities which would seem to destroy the viability of these theories. However, 



singularities are not necessarily disasters. Some of them may be hidden behind horizons [[16 
Alternatively, they may be resolved by stringy considerations. 



Recently, Polchinski and Strassler considered another model [H], the M = 1* model, 



in which mass terms for the three chiral multiplets that make up an M = 4 multiplet are 
included. Several interesting generalizations of |17]] have also appeared fT8H23|. At high 



energies, the theory is the usual M = 4 model, while at low energies, the supersymmetry is 

broken to M = 1. Unlike gauge theories with higher supersymmetries, M = 1 theories can 

be confining. In fact, this gauge theory has a large class of degenerate but discrete vacua 

2TI.EH] which include Higgs, confining, and oblique confining phases. In these cases there 



is a mass gap, but there are also vacua with unbroken U(l) subgroups and thus massless 
photons. 

Polchinski and Strassler constructed the supergravity dual of this model as follows. They 
assumed that a 3-form field strength is turned on whose strength is proportional to the mass 
m of the hypermultiplets. The field strength induces a Myers dielectric effect on the D3 



branes, essentially expanding them out in the transverse directions to a two sphere |26| , |27 
The two sphere is effectively a D5 brane, or one of its S-duals, with n units of magnetic 
flux. These flux units minimally couple to the 4-form Ramond-Ramond gauge potential, and 
hence they correspond to D3 brane charges. The background metric, for large distances, 
that is for the UV limit, is dominated by the D3 brane charges. But as one moves toward 
the IR, one starts probing near the D5 brane and here, the metric is dominated by the D5 



brane. In between the metric interpolates between these two limits. 

The presence of this D5 brane significantly changes the physics |T7|]. For example, at 
large separations, magnetic charges are effectively confined, since a D string will have finite 
tension as it approaches the D5 brane. However, the electric charges are screened, since 
the fundamental string can end on the D5 brane. Likewise, for the S-dual picture, electric 
charges are confined since a fundamental string has finite tension near an NS5-brane, but 
the magnetic charges are screened because the D-string can end on the NS5-brane. 

Polchinski and Strassler also demonstrated how the supergravity duals map into the 
various degenerate vacua of the M = 1* theory. The D3 branes can expand into more than 
one 5 brane and in the regime where the supergravity limit is valid, one can find a one to 
one correspondence between the various ways of dividing up D3 branes into 5 branes and 
the degenerate vacua. For these vacua, one expects a gap, in other words, the absence of 
massless particles. Actually, this is not completely true, since many vacua have unbroken 
U(l) gauge groups. However, the number of massless states is at most of order \rN and 
their effects are not seen in supergravity duals. 

If one were to consider the M = 1* theory at finite temperature, then the multiple 
vacua and gaps have profound effects. Consider first the entropy. The entropy of a massless 
theory in four dimensions scales as T 3 . However, a theory with a gap will have its entropy 
suppressed exponentially when the temperature falls below the mass gap scale. In the large 
N limit, this should occur as a phase transition between high and low temperature phases, 
where at some critical temperature the free energies of the two phases is equalized. Another 
effect of finite temperature is that the degeneracy between the separate vacua should be 
lifted, leading to many metastable states. As the temperature is raised above the mass gap 
scale, these metastable vacua should disappear altogether, and at very high temperatures, 
the theory should behave just like the M = 4 theory at finite temperature. 

In the supergravity dual of these theories, the high temperature phase is described by 
a black hole in an asymptotic AdS space. At high temperatures there is a unique vacuum. 
The multiple vacuum configurations of 5-branes described by Polchinski and Strassler will 
now have a higher free energy than the vacuum. The 5-branes were introduced to account 



for naked singularities that inevitably appear and which were specifically discussed in p8 
However, at high temperatures, singularities are shielded by the horizon; hence, there is no 
need for the 5-branes. 

For the high temperature phase, if one considers a spherical 5 brane probe with some D3 
brane charge, then one should find that the free energy is minimized when the 5 brane lies on 
the horizon. One should also find that the energy is the same as if there were only D3 branes 
and no 5 branes, since the 5 brane has no net charge, and hence would be undetectable at 



the horizon. 

However, in the high temperature phase, there still may exist local minima for the 5- 
branes outside the horizon. These configurations correspond to metastable vacua, valid in 
the probe limit described in [H], and which can be in a partial Higgs, confining or oblique 



confining phase. As the temperature is increased, these local minima eventually disappear. 

A crucial step for going beyond the probe calculation is the realization that the probe 

calculation itself is unchanged for an arbitrary distribution of parallel D3 branes [|17|. Hence, 



one can construct a spherical shell of D3 branes, of order N in number, and still use the probe 
calculation to determine the radius of the shell. This is not to say that the exact supergravity 
solution has been constructed, since one still needs to compute the back reaction. But it 
does allow one to build an interpolating metric between a D3 brane geometry and a 5-brane 
geometry. However, in the finite temperature case, one cannot use this approach as a starting 
point, since a configuration of D3 branes outside the horizon is not a solution to Einstein's 
equations. 

One can also explore the effects of the hypermultiplet mass on the entropy in the high 
temperature phase. If the hypermultiplets are given a small mass m, then there are polyno- 
mial corrections to the entropy. In particular, to order m 2 the free field limit of the entropy 
is given by 

S = —N 2 T 3 - -N 2 m 2 T + 0(m 4 ) (1) 

3 4 v ; v ; 



The strong coupling calculation for the massless case was carried out in [^J , where the area 
of a blackhole in AdS space was computed. These authors found that 

S = ^N 2 T\ (2) 

a factor of 3/4 from the free theory result. An immediate question is to find the massive 
corrections in the supergravity limit and compare this with the free field result. 

In this paper we study the finite temperature effects outlined above for the supergravity 
dual of the J\f = 1* theory. Our starting point is the black D3 brane solution |3(| of Type IIB 



supergravity perturbed by a 3-form gauge potential. A 5-brane probe calculation quickly 



reveals several features of finite temperature. First, the T = minima Jl7j for 5-branes 
wrapping 2-spheres of particular nonzero radii are now only local minima; the true minimum 
appears at the horizon, r = r^ where the free energy is minimized. The local minima, which 
can be thought of as metastable vacua, move to smaller radii and eventually disappear as 
the temperature is increased , but at temperatures that can be much higher than when the 
field theory is weakly coupled. We demonstrate that at the horizon a wrapped 5-brane, 
which carries n units of D3-brane charge, cannot be distinguished from a configuration of n 
D3 branes. 



Although the major thermal effects in the probe calculation can be obtained quite simply, 
other questions require the finite temperature modification of the fields of the Type IIB 
theory through order m 2 . We therefore compute exactly the 3- form field strength in the 
black D3 brane background. The metric, 5-form field strength and dilaton are then sourced 
at order m 2 by bilinears in the 3-form, and we solve the equations which determine the 
back reaction to this order. Given the modified metric, the relation between temperature T 
and horizon distance r# is recomputed, and we find that the horizon shrinks for fixed T by 
an amount of order Ar# ~ m 2 R 2 /T where R is the AdS scale. We presume that a more 
complete calculation would show that Th — > at some temperature T c which is then the 
critical temperature for the phase transition in the theory. 

The computations for the back reaction can be specialized to the T=0 case, where 
we must require that the homogeneous modes agree in their 5*0(6) representation and r- 



dependence at the boundary with fields of the dimensionally reduced IIB supergravity |]3T 
In general the order m 2 inhomogeneous solutions approach the boundary at different rates 
than the homogeneous modes. For example, the dilaton perturbation equation is that of the 
second Kaluza-Klein excitation, the A = 6 field in the 20' representation with an order m 2 
source. The 20' representation agrees with the general relation between the D = 10 and 



the D = 5 dilaton postulated in p3[, but the solution does not agree with the order m 



term in the interpolating dilaton of [ I7| . The back reaction equations of metric components 
transverse to the D3-branes describe scalars with A = 8 in the singlet representation and 
A = 6 in the 84 representation of 5*0(6). We also show that the back reaction on the metric 
and 5-form lead precisely to the order m 2 terms in the probe action, which were inserted in 
j jTTfl as a requirement of supersymmetry. All of these T = results are a useful check that 
the rather complicated back reaction equations are correct. 

We also compute the decreased area of the horizon in the (approximate) modified geom- 
etry and thus obtain the entropy 
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S = —N 2 T S - 0.17UN 2 m 2 T + 0(m 4 ) (3) 

including the order m 2 correction. Comparing with (HD we see that the large N correction 
is 23% of the free field result. 

Technically speaking, it is noteworthy that finite temperature effects can be computed 
essentially analytically as far as we have gone, and this raises the hope of going beyond 
the leading order m 2 approximation. At T = this proved to be easy W7\ since the probe 



calculation is not altered when the D3 branes are expanded around 2-spheres. This fortuitous 
simplification is unlikely at finite T, since the only solution to the equations of motion when 
m = has all D3 branes at the origin. 



In Section II we consider the probe calculation. In Section III we compute the linear 
solutions for the 3-form field strengths in a finite temperature D3 background. In Section IV 
we obtain the back reaction to the dilaton field. Order m 2 temperature dependent corrections 
to the metric and 5-form are discussed in Section V, and the entropy calculation is presented 
in Section VI. The appendix contains a more detailed calculation for the disappearance of 
metastable vacua. 

II. PROBE CALCULATION FOR THE VACUA 

In this sectionF] we look for nontrivial solutions for a D5 brane probe in a finite tempera- 
ture D3 brane background geometry. Temperature is introduced via the thermal D3-brane 
solution |3(J of Type IIB supergravity. The metric is 

ds 2 = (Z(r))-^[-f(r)dt 2 + dx 2 ] + Z(r)*[-^r + r 2 dn 2 5 ] (4) 

fin 

Z(r) = £ (5) 

m = i - ^ (6) 

where R 4 = 4ngNa' 2 . The solution has a horizon at r = r#, which corresponds to the 
temperature T = r#/(7ri? 2 ). The accompanying dilaton, axion and self-dual 5-form are 
given by |T7 



e- = 9 (7) 

C -i <"> 

F 5 = d X4 + *dX4 (9) 

dx° A dx 1 A dx 2 A dx 3 , , 

* = IW) — • (10) 

While thermal effects cancel out for the * operation contained in F 5 , they will play an 
important role when * acts on other field strengths. 

The action for the D5 brane is the sum of Born-Infeld and Chern-Simons parts, and is 
given by 

1/2 



S = -^- f d 6 t[-det(G ll )det(g- 1/2 e* /2 G ± + 27ra'F)] + /i 5 f(C 6 + 2ira'F 2 A C 4 ) , 
g J J 

(11) 



x the notation of [[D]] is used throughout the present paper 



where 

2itaT 2 = 2na'F 2 - B 2 . (12) 



As in |T7||, G\\ refers to the pullback of the metric along the R 4 , while G± is the pullback of 
the metric onto the S 2 on which the 5-brane is wrapped. C^ is the Ramond-Ramond gauge 
potential of the D3 branes, while C@ and B 2 are potentials for the 3-form perturbation of 
the background dual to the mass operator m tr($f + $2 + ^3) °f the M = 1* field theory. 

The field F 2 is the U(l) gauge field strength of the wrapped 5-brane, and is assumed to 
contain n units of flux, so that 

f F 2 = 2im. (13) 

Since F 2 couples to C^ in the Chern-Simons term of (|TTD, n is the effective D3 charge of the 
5-brane. As in [pH , we assume that the first term in (|TTD is dominated by F 2 . The condition 



n 2 » gN must then hold in order for the probe calculation to be valid ]T7 . 

In the zero temperature limit, the leading order n/y/gN term from the Born-Infeld part 
of ( |ITD cancels with the C^ term in the Chern-Simons part. The radius of the S 2 on which 
the D5 brane resides is then chosen to minimize the sum of first order correction terms. 
At finite temperature, leading terms no longer cancel because det G» is modified by the 
Schwarzschild factor |[32| , viz. 

det(G||) = Z- 2 (l-rVr 4 ). (14) 

There is then a residual r-dependent term which must be added to the non-leading terms 



computed in |T7] . The sum of all these terms is then minimized to find the favored position 
of the probe. The modification of detG\\ is the major thermal effect on the probe calculation. 
Thermal modification of other quantities is considered in later sections of this paper and 
produces corrections of higher order in r 2 H /r 2 which can be ignored when the 5-brane is far 
from the horizon. The near-horizon effect of these corrections will be included qualitatively 
here and then supported by the work of later sections. We now proceed to implement the 
probe calculation just outlined. 



At finite temperature, the contribution to the action from the Born-Infeld term in (11) 
which is of leading order in the expansion parameter n/y/gN is 

27 f^/^f-dct( G|l) V /2 (det(^^/ 2 



5o = -^^/^(-det(G l |)) 1/a (det(F J ,)) 1 



= -^w^vyr^^M (15) 

where r# is the position of the horizon and V is the volume of R 4 . The leading contribution 
from the Chern-Simons term is 



2 W/i 5 />, A C, = +^4vr V^V. (16) 

J g W> 

Dominant thermal corrections are absent since the metric does not appear in this term. 
Further only the parallel components of the Ramond-Ramond gauge potential contribute 
to the integral, so we have set C 4 = X4- The zero temperature first order correction to 
the Born-Infeld term plus the contribution of the Cq plus the order m 2 term required by 
supersymmetry is given in [O] as 



AS = -^r 2 (r-r )V, (17) 

ga'n 

where 

r = ir\m\na'. (18) 

Adding the terms in (|15D , fllBD and ([H]) leads to the expression 

-v=Mm r,{ J TZ ^- i)+r2(r - r ° )2 )- (i9) 

This result for the free energy of the probe is valid if r >> r#. Note that for r# = 
(|19|) has a minimum at r = r and a maximum at r = r /2. We now discuss thermal effects 



in this region of r. For rn « r the minimum and maximum shift by 



im r] 
2gWrl 



7m r H 



327m 2 r 8 H 
gNrb 

Therefore, as th is increased, the maximum and minimum are pushed toward each other, 
and the minimum disappears when Ar max ~ r . This occurs at the critical horizon size 
r H ~ (gN/n 2 y/ 8 rQ which is well inside of r and thus within the region of validity of the 
calculation. 

To understand what happens closer to the horizon we must incorporate thermal effects 
on the second term in ([19]). One effect is that the r 4 term in fll7]) should be multiplied by 



1 — rfj/r 4 since it comes from the Born-Infeld part of (|Tl|)- There are other effects from 
the 3-form perturbation and its back reaction on the metric and 5-form. The linear term in 
r comes from the Chern- Simons term 

Ms f(C 6 -B 2 AC 4 ), (22) 



and goes to zero at the horizon due to a conspiracy between Cg and B%. The r\ term 
has contributions from both the Born-Infeld and Chern-Simons parts of the probe action. 



The Born-Infeld contribution is again suppressed by the factor Jl — rjj/r 41 , but the Chern- 
Simons term also vanishes since the order m 2 corrections to C 4 approach zero. These 
statements will be verified to order m 2 in sections 3 and 4. 

The conclusion is therefore that the 5-brane has lower free energy at the horizon, so that 
the minimum near r = r$ is a metastable state even for small finite th and is a degenerate 
vacuum state only if r# = 0. There is still another thermal effect which we discuss in section 
6. Due to the back reaction of the 3- form perturbation on the metric the horizon size itself 
shrinks as a function of m 2 for fixed black hole temperature. Presumably Tb vanishes at 
some finite temperature T of order T ~ m. This value would then be a critical temperature 
for the probe calculation. In the low temperature phase T < T Q the 5-brane is stabilized 
near r = ro; in the high temperature phase T > To it retreats to the horizon. Work is 
continuing to verify this suggested picture. 

Standard no-hair theorems should imply that the D5 brane at the horizon is equivalent to 
n D3 branes, since there is no net D5 brane charge. At the horizon, there is no contribution 
from the Born-Infeld action. Hence, one would need to verify that the actions coming 
from the Chern-Simons pieces are equivalent. The difference between the two actions is 
proportional to (p2|). As previously stated, the order m 2 contribution will be shown to be 
zero at the horizon. 

The generalization of the previous analysis to NS5 probe branes with p units of D3 brane 
charge is straightforward. To proceed, one describes the D3 background using the S'-dual 
description 

, ar + b , . 

cr + d v ' 

with field transformations 

g' = g\M\ 2 G' MN = G MN \M\ C 4 = C 4 

G' 3 = G 3 B' 6 - t'C 6 = (B 6 - rC 6 )M- 1 (24) 

where M = ct + d. Hence, the leading order term in ([19|) is invariant under this transfor- 
mation. The first order terms have the same transformation as in 11711, 



AS = ^r 2 (r-r ) 2 V. 
g 6 a'p 

r = ng\m\pa' (25) 

Hence, the total action that is minimized is 



S 2/i 5 firgp 



2 



-/" 



'1/ 



l-r^/r 4 -l)+r 2 (r-f ) 2 • (26) 



V g 3 a'p \ 2N 

Finally we point out that a more precise analysis in the Appendix shows that the 
metastable D5 probe solution near r = r^ disappears above the temperature 

/ n 2 \ 3/8 

T^k[ \m\ (27) 

\47cgNj ' ' V ; 

1 / ^\A \ /® "I /H 

« = — ( — ) (3351797+171(2659)^) ^0.5522 
16 \7rV v ' 

Likewise, for NS5 brane probes one finds that the solution disappears at the temperature 



T^ K (^-) \m\, (28) 

\4nNJ ' '' V ; 

Hence these metastable probe minima can survive well above the critical temperature, 
at least for those probes that satisfy the probe condition n 2 » gN. This contrasts with 
the weak coupling analysis. At weak coupling, the presence of temperature introduces the 
effective mass term T 2 tr(<f) 2 ) to the effective potential. Hence, there are metastable Higgs 
vacua here as well. However, in this case these local minima are washed out when T ~ m. 

What is intriguing about these metastable vacua is that one could have solutions where 
only some of the D3 brane charge is outside the horizon. The corresponding vacuum is a 
metastable state where say part of the SU(N) gauge group has been Higgsed, but the rest is 
unbroken and unconfined. Likewise, we can have metastable vacua where part of the gauge 
group is confined, but the rest is unbroken and unconfined. 

III. THE 3-FORM PERTURBATION 

As in |T7| the first step in the study of thermal effects is to obtain linear perturbations 



of Hg and F$ = Fg — CHg which are dual to fermion mass terms in the boundary gauge 
theory. Specifically we must solve the linearized equations of motion and Bianchi identities 
in the background (@), (0)-©: 

(29) 
(30) 

(31) 



(32) 
(33) 









d* F 3 


= F 5 AH 3 






d * (-H 3 
9 


-gCF s ) 

dF 3 


= -gF 5 A F 
= = dH 3 


ing the fields 


as 


the complex 


combinations 






G 


3 — F 3 — 

T = C + 


tH 3 
i/9, 



10 



the equations of motion in ( p9"D can be recast into the compact form [17| 

d * G 3 + igG 3 A F 5 = (34) 

dG 3 = (35) 

The black hole metric affects only the radial dependence of G 3 , the dependence on angles 
must be that of the lowest spherical harmonic on S 5 in order to be a source for the fermion 



mass in the 10 or 6/10 representation of SO (6). Therefore we postulate that G3 is the 
combination 

G 3 = a(r)T 3 + (3(r)V 3 (36) 



of the same 3- forms used in |17|| . The 3-forms T 3 and V 3 are constructed from a constant 
antisymmetric tensor T mnp which is either self-dual or anti-self-dual corresponding to the 10 
or 10 representations, respectively. These forms are simply written in terms of Cartesian 
coordinates y m , m = 1, 2, • • • , 6 in the space perpendicular to the 3-branes, with r 2 = y m y m ; 

T 3 = ^T mnp dy m A dy n A dy p 

S 2 = \T mnp y m dy n A &tf (37) 

V 3 = d(\nr) AS 2 

One can easily derive the properties 

*6 T 3 = = ±iT 3 
* 6 V 3 = ±i{T 3 - V s ) 

dT 3 = (38) 

do 2 == 3-Z 3 

dV s = -3d(lnr) AT 3 

where *6 is the Poincare duality operation in flat R 6 . 

We now proceed to solve the equations (|34]) . Substitution of the ansatz (|36D in the 
Bianchi identity readily gives 

dG 3 = (a'(r) - 3p(r)/r)dr A T 3 = (39) 

which determines the relation /3(r) = ra'(r)/3. One can then directly show that 

G 3 = l -d{a(r)S 2 ). (40) 



11 



The equation of motion in fl34|) contains the duality operation with respect to the full metric 

cam, 



JL qq'ss'tl/ 

i c mnpqsty y y 

where y mn is the transverse piece of the metric of (H) rewritten in terms of y 



* Gs = ^r e mn P qstg qq 'g ss 'g tt 'G q/s/t/ dy m A dy™ A dy? A ^(r) X4 (41) 



III 
1 



f, y m y n \ f -i, ,y m y n ,, , 

When compared to the case of zero temperature (r# = 0), one sees that the radial compo- 
nents of G 3 dualize with an extra Schwarzschild factor f(r) coming from g rr . 

To handle this simply we split the forms T 3 and V 3 into radial and angular parts, using 

dy m = y m d(lnr) + u m (43) 

where io m is a set of angular 1-forms satisfying 

du m = -dy m A d(lnr) (44) 

Then we have 

T 3 = f 3 + f 3 (45) 

with 

f 3 = l -T mnp y m d(\nr) Au n Au? 

= ^T mnp y m d(\nr)Ady n Ady p 

= V 3 (46) 

f 3 = ^T mnp u m Au n Au p 

Note that V 3 = <i(lnr) A Sg = V 3 has no angular part. The differentials are 

df 3 = -3d(\nr)AT 3 df 3 = dV 3 (47) 

Under flat * 6 duality, we have 

* 6 (f 3 +f 3 ) = ±i(f 3 +f 3 ) (48) 

so that *qT 3 = ±iT 3 and * 6 T? = ±iT 3 . Under the duality ( f4"T| ) of the 10-dimensional metric 
(El), one thus has 



12 



*f 3 = gf(r) * 6 f 3 A X4 = ±igf(r)f 3 A X4 

*f 3 = g * 6 f 3 A X4 = ±WT 3 A X4 

For V 3 (|38|) implies 

* V 3 = * f 3 = ±igf(r)f 3 A X4 

Combining these duality relations with the ansatz (^), we find 

* G 3 = a{r) * (f s + f 3 ) + (3{r) * % 

= ±i[a(r)(f(r)f 3 + f 3 ) + f3(r)f(r)f 3 ] A X 4 



(49) 



(50) 



Finally we note that 

G 3 AF 5 = G 3 A(d X4 +*dX4) 



d{G s Ax4. 



(51) 



(52) 



where we have used the fact that G 3 A *dx4 = and dG 3 =0. 
Thus the equation of motion can be simply rewritten as 



d[*G 3 - %G 3 A X4 ]=id(± [a(r)(f(r)T 3 + T s ) + (3(r)f(r)T 3 - (a(r)T 3 + f3(r)T 3 )\ A X 4 , 
= (53) 



Using 



d X4 = (Z'/Z)dr A X4 = -4d(lnr) A X4 



(54) 



and (|47D it is straightforward to work out the differential, substitute (3 = ra'/3, and obtain 
the differential equation for a(r): 



rf(r)a"(r) + (ll/(r) + rf(r))a'(r) + 3(/'(r) + 7/(r)/r + 7r)a(r) 







(55) 



where 7 = —7 in the self dual case *qT 3 = %T 3 and 7 = 1 in the anti-self-dual case 
*qT 3 = —iT 3 . We now restrict to the latter case for which the fluctuation G 3 is in the 
10 representation dual to the fermion mass term trXX in the M = 1* field theory. Inserting 
the specific form of /(r) and making the substitutions a(r) = F(w)/r 4 with u = rjj/r 2 , one 
finds the Legendre equation 



F"(u) 



2u 



1 — u 



;F'(u) 



1 



4(1 - u 2 



-F{u) = 0. 



The general solution is 



a{r) = -j[aPi(r 2 H /r 2 ) + bQ_ h (r 2 H /r 2 )} 



(56) 



(57) 
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Only the first term is regular at the horizon, so we set b = and observe that this solution 
has the expected JT7J behavior 1/r 4 as r — > oo associated with the fermion bilinear trXX 
which is an operator of scale dimension 3. The coefficient a is set by comparing ([57]) to the 
zero temperature result |T7J, giving the solution 



a(r) = -^-L—P_ h (r> H y). (58) 

We also observe that there is no regular solution which approaches the boundary at the rate 
1/r 6 associated with a vacuum expectation value (trAA). 

In the zero temperature case, G 3 is singular at the origin. It was thus necessary to insert 
D5 or NS5 branes wrapping two-spheres in order to interpret the singularity. In this case, 
where there is a black hole with a finite Schwarzschild radius, all singularities are hidden 
behind the horizon. Hence, one is no longer compelled to insert 5 branes. From the gauge 
theory point of view, one should interpret this as being in the high temperature phase where 
the system acts like a finite temperature M = 4 theory. Only when r# shrinks to zero will 
the solution for a(r) be singular, and thus necessitate the insertion of 5 branes. 

Let us now turn to the 6-form dual potentials. Making use of the Bianchi identity and 
the relation of V 3 to T 3 , we can rewrite the expression for G 3 in (|3"6j) as 



G 3 = a(r)T 3 + (a(r) + (l/3)ra\r))f 3 . (59) 

The equation (|53|) implies that *G 3 — iG 3 AX4 is a closed 7-form. Hence, it can be written 



as 
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* G 3 - iG 3 A X4 = ~d{B 6 - rC 6 ) 

% 



9 -d f^ ° (r) + /(r)( ^ r) + ra ' (r)/3) S 2 ) dx° A d X ^ A ^ 2 A dx\ (60) 



where the last line is derived using the relations in fl38|) . If we assume that C — 0, and hence 
t is imaginary, then Cq is given by 

C 6 = -| (a(r) + /(r)(a(r) + ra'(r)/3)) lvaS 2 A X4 . (61) 



It is instructive to compare this result for Cq to the NS-NS 2- form. Using ( p2|) and (fJ0|), 
one finds that 

B 2 = -|a(r)Im^. (62) 

Since /(r) = at the horizon, we learn that C$ = B% A \4 at r = th- This confirms the 
claim in the previous section that a D5 brane with D3 brane charge n has the same action 
at the horizon as n D3 branes. 
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IV. BACK REACTION OF THE PERTURBATION ON THE DILATON 

In this section we compute the order m 2 corrections to the dilaton due to the temper- 
ature dependent G3 perturbation of the previous section. The modified dilaton actually 
contributes only higher order corrections to many of the physical effects we are interested 
in, but the calculation is simpler than for the 5-form and metric and serves as a prototype. 
The result is of interest for T = 0, since the order m? correction to the dilaton, and hence 
the gauge coupling, is independent of the vacuum state. 

We must solve the equation 

D M D M * = -±H MNP H MNP + ^F MNP F MNP , (63) 

which was obtained from the more general dilaton equation of motion in |17]] by specializing 
dilaton and axion sources to their constant values in (Q) and @. The metric is that of (|J) 
with perpendicular part rewritten in terms of "Cartesian" coordinates y m as in (f42|). The 
source can be expressed in terms of G3 as 

D M D M ® = ^g mm 'g nn 'g pp '{G mnp G m , n , p , + c.c] 

2AZ2 r 2 

in which lower pairs of indices are summed over their 6 values. Using the representation 
(|36|) and the fact that y m V mnp = y m T mnp , we obtain 

/"f fi I zSL \ i Lii m n ii m ' n 

f ^ T mnp^ Jr mnp "T 2 " ^ J mnpy ^ J m'np 

Ot y ') J- pqr J- pqr 5 V V ^-mpq^-npqi \F* I 

where a(r) is the same functions derived in the previous section. 

We now need the specific form of the anti-self-dual tensor T mnp whose non-vanishing 
complex components are fTT| 



-'■pqr -'■pqr -'pqr f^^-pqr \~®) 

where e pqr is the Levi-Civita symbol in 3 (complex) dimensions. By conversion to complex 
coordinates and reconversion, we obtain 

T T =0 

- 1 pqr - 1 pqr w 

y m y n T mpq T npq = 2m 2 r 2 Y(y l /r) (67) 

(y 1 ) 2 + (y 2 ) 2 + (y 3 ) 2 - (y 4 ) 2 - (y 5 ) 2 - (y 6 ) 2 



Y(y*/r) 
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The function Y(y l /r) is a harmonic polynomial for the 20' representation of 50(6) and has 
"eigenvalue" — 12/r 2 of the flat 6- dimensional Laplacian in the y l . Thus the final equation 
for the lowest order dilaton perturbation is 

777 CI 

D M D M § = --4[« 2 (r) - /(r)(a(r) + ra'(r)/3) W/H (68) 

where we have assumed that m? is real to prevent a source term for the Ramond-Ramond 
scalar field. 

The next step is to find the specific form of the Laplacian. Since the source depends only 
on r and angles, we can drop t and x derivatives. The solution must also be proportional 
to the harmonic Y(y l /r) so that angular derivatives just give the eigenvalue —12. Thus 

D M D M ® = D M D M 4>{r)Y{y i /r) = ^V/dr 2 + (- + ^)d/dr) - -^}<P(r)Y (y l / r) 

Mr r j r z j[r) 

(69) 



In the limit of zero temperature (/(r) — > 1) (|69| ) and (38) combine to give the radial 
differential equation 

This equation has homogeneous solutions <f>(r) oc r 2 , r~ 6 which are the power laws of the 
irregular and regular solutions for a field of scale dimension A = 6 in the 20' representa- 



tion of 5*0(6). This is the second Kaluza-Klein excitation of the dilaton/axion |31 ]. The 
inhomogeneous solution is 

,. . m 2 R 4 ,_,. 

^ r = ^7^' 71 

2r z 

Presumably this nontrivial behavior for the dilaton describes the scale dependence of 
the field theory gauge coupling due to the introduction of mass for the chiral multiplets. 
However, we do not clearly understand the angular dependence of $ = <f)(r)Y {y l /r) and thus 
cannot give a more precise interpretation. For a softly broken M = 4 to M = 2 theory one 
finds angular dependence for the effective couplings on the coulomb branch [Q. However, 
for M = 4 broken to M = 1 there is no coulomb branch, so the Wilsonian interpretation 
seems less clear. 

In this connection it is interesting to note that for the (T = 0) solution a(r) ~ 1/r 6 , 
which is dual to the field theory vev (tr(AA)), the source term in (RSI) vanishes. Thus the 



dilaton does not run in a situation where no masses are turned on [041133 



In W7\ an interpolating function for the dilaton was given that followed a 5-brane geome- 



try in the IR to a 3-brane geometry in the UV. For example, for a vacuum that corresponds 
to one D5 brane with D3 brane charge N, the interpolating dilaton was given by 
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2$ 2 

e =g 



P- 



Pt + Pi 



(72) 



where 



P- 

Pc 



Vi + 2/2 2 + y 2 + 



R 2 



Vi 2 + V5 2 + V6 2 - r 
r = Tia'mN. 



If we expand this for large r we find that the order m 2 dilaton correction is 

k^ ( ^ m2Ri 



(73) 



(74) 



This is an 50(6) singlet, hence the first order correction in the interpolating dilaton is not 
quite consistent with the supergravity equations of motion. If one were to have considered 



an NS5 brane, then one would have found ( [74]) but with the opposite sign. In fact, if one 
were to solve the supergravity equations but with 5-brane sources inserted, one would find 
that their effects are of higher order in m. The order m 2 term should be independent of the 
particular vacuum the M = 1* theory is in. 

At finite temperature the equation in (|68|) is more complicated but tractable. Using the 
variable z = Z(r) = R A /r A and setting k = rjj/R 4 , we find the radial differential operator 



[d 2 /dr 2 + (- + ^Ad/dr - 
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16- 



212 



d 2 



k d 



r f J r 2 f(r) * R 2 dz 2 1 — kz dz 4z 2 (l — kz)' 



which is essentially hypergeometric. Its zero modes are: 

3 3 3 3 



Mz 






7lk 2 Z2 



whose Wronskian is 



w(4>w 



s 



nk 2 (l 



» (75) 

(76) 
(77) 

(78) 



The limiting form of the solutions at the boundary (z -^ 0) is also given in ([76]). At the 
horizon (z — > 1/k), the function <p\{z) has a logarithmic singularity while <f>2{z) is regular. 



The finite temperature equation with source fl68 ) can be rewritten as 

k d 3 S(z) 



^dz 2 

S(z) 



1-kzdz Az 2 {l 
m 2 g 2 R 2 



32z 7 2 



kz) x 1 — kz 

a 2 (z) - f(z)(a(z) - Aa'(z)/3) 2 



(79) 
(80) 
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The method of variation of parameters gives the solution 



rl/k rz 

0i (z) \ dz<p 2 (z)S(z) +(p 2 (z) dz(f) 1 (z)S(z) 
Jz Jo 



(81) 



which is regular on the horizon and vanishes on the boundary at the same rate as ( |71~D . This 
solution does not appear to be integrable in closed form. 



V. BACK REACTION TO THE METRIC AND SELF DUAL 5-FORM 

In this section we compute leading order corrections to the metric and self dual 5-form 
when the field G3 is turned on. The Einstein equations are 

1 q 2 q 2 q 2 

Rmn — -^QmnR = —^FmpqrsFn + -tt(GmpqG* n + c.c) — —gMNGpQRG* , (82) 

where the contribution of the dilaton is ignored since it affects only higher order terms. 

Let us consider the G 3 terms first. Using the results in fl32|), ( f4T3| ) and (|42|), it is straight- 
forward to show that 



|2 / /„,m„,n 






24Z 



(a(r) 2 (3I mn - 2W mn ) + (a{r) + ra'(r)/3) 2 (^-f- + 2f(r)(W mn + J mn ))) , (83) 



where a{r) is the same as in (ESF), |T 5 | 2 = T pqr T* pqr and 



pqr 

m ,n 



1 / y m y 

- d. 



1 mn .., \ - 1 1 1 1 > 9 



5 V T A 



3 f m m * y k y l 

Wmn = p 

Using the fact that 



W mn = -— lT mpk T* npl ^- + c.c\- I mn . (84) 



y{T mp qT* pq + c.c) — -5 mn \T 3 \ , (85) 

one finds that W mn is traceless. We also see that W mn is strictly an angular tensor since 
y m W mn = because of the antisymmetry of T mnp . 
Taking the trace of the expression in (p3|) gives 



G mpq G* m ™ = i^ (a(r) 2 + (a(r) + m'(r)/3) 2 ) . (86) 

Hence the contribution to the energy-momentum tensor from the G3 kinetic terms is 



q 2 Q 2 

n 2 \To\ 2 / ?/ m ?/ n \ 

^|^(a(r) 2 - /(r)(a(r) + m'(r)/3) 2 ) (/ mn - 4W mn - r\r) V -^-) . (87) 

for the orthogonal components and 

- ± 9ij G pqr G*^ = -M («(r) 2 + (a(r) + ra'(r)/3) 2 ) % (88) 

~9uG pqr Cr"' = M ( a(r f + ( a ( r ) + ra '(r)/3) 2 ) /(r), (89) 

for the space-time components. 

We next compute the corrections to F 5 = F 5 — |(C g A H 3 — B 2 A F 3 ), which in terms 
of Gs is 

F 5 = F 5 + |(</> A G* 3 - r A G? s ), (90) 

where t/> = ^a(r)S 2 - The equation for motion for F 5 is 

F 5 = *F 5 . (91) 

There are two correction terms that need to be added to the lowest order term 

F° 5 =dx 4 ,o + *dx4fi (92) 

in order to satisfy (|9lD , where \4fi is the 4-form field in (|i~0l). First, since there are corrections 
to the metric, Xi,o is no longer harmonic and so the expression in ( |92"1) is not closed. We 
will return to this point once we discuss corrections to the metric. Second, the improvement 
term in ( |90"D is not self-dual, so it is necessary to add another correction term to F 5 to 
compensate for this term. A good guess is to add the dual of the improvement term to 
Fg, which would make the sum of this term and the improvement term manifestly self-dual. 
However, one first needs to verify that the dual is closed. To see this, we note that 

G 3 = ^a'(r)drAS 2 + a(r)T 3 

= (-a'(r)dr + r~ 1 a(r)') A S 2 + -r 3 a(r)d(r~ 3 S 2 ), (93) 

where the second term on the rhs of the second line in ( p3[ ) has no radial components. 
Therefore, we find that 

|(V A G* 3 - r A G 3 ) = ^a 2 (r) (S 2 A T* 3 - S% A T 3 ) , (94) 



with no radial components. Therefore, the dual of (|94l ) is 
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" /,r( ^ * 6 (S 2 A T* - S* 2 A T 5 ) A dx° A dx 1 A dx 2 A dx 3 



12Z 2 

<?a(r) 2 |T 5 | 2 



rdr A dx° A dx 1 A dx 2 A <ix 3 , (95) 



72Z 2 

which only has a radial component in the orthogonal directions. Since a(r) and Z only have 
r dependence, the expression in (|95|) is clearly closed. 

There is an m 2 correction coming from the cross term of the improvement term and its 
dual with F®. For the angular, i. e. magnetic, components, the contribution of this cross 
term to fl32|) is 

ig 2 a 2 {r)rZ' 6 vVfofrp T * \_o(t T* W- 

288Z 2 €i i kml( l r 2 \ \ Pi k nlq C - C -J L \ l VJri^klq C.C.J j — 

9 2 a 2 (r)\T 3 \ 2 f y m y n \ 

36Z - ^--H-'J (96) 

where we have used the fact that *§T 3 = —iT 3 . It then follows that the electric cross term 
piece is 

9 2 a 2 (r)\T 3 \ 2 

36Z 2 9MN > (97) 

where gMN is the AdS$ Schwarzschild metric given in ([|) 

We can now combine the results of (|87|), (|S8|), (|S9]), (|9~6|) and (|9T| ) into one source term, 
Jmn, for the right-hand side of (^). For the orthogonal components we have 

„2|T^ 12 / ii rn ii n \ 

Jmn = ~ 9 -^ (« 2 (r) + 3/(r)[a(r) + m'(r)/3] 2 ) (hl mn - f-\r) V -^- J 



21/71 |2 



^|T ; 



12Z 
while the space-time components are given by 



(a 2 (r) - /(r)[a(r) + ra'(r)/3] 2 ) (W mn + I mn ), (98) 



J<, = + 9 -^- (a 2 (r) - 3/(r)[a(r) + m'(r)/3] 2 ) % - (99) 



<? 2 |?V 2 



J °° = -UIT- (" 2(r) ^ 3 ^ r )H0 + m'(r)/3] 2 ) /(r). (100) 



Upon inspection of fl9~8"[) one sees that the right hand side consists of three types of terms. 
The radial term y m y n jr 2 and the two angular terms I mn and W mn . Hence, the corrections 
to the transverse metric components should be of this same type. Likewise, (^9j) indicates 
that the spatial piece of the metric can have a correction proportional to rjij while (|100|) 
indicates that the temporal component has a correction proportional to f(r). Hence, we can 
write down the following ansatz for the metric 

ds 2 = - [Z- 1 ' 2 + ho(r)) f{r)dt 2 + (Z~ l l 2 + h t (r)) r ]ij dx i dx j + 

1 dy m dy n . (101) 



7 / II 






(5Z 1 / 2 + p{r))I mn + r\r){Z 1 ' 2 + q(r))^- + w{r)W n 
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Let us now return to the corrections to Fg induced by metric corrections. Let us write 
the 4-form x^( r ) = X4,o( r ) + X4,i( r ), where X4,i( r ) is assumed to have components in the 
four space-time directions, 

X 4 ,i(r) = xi(r)dx° A dx 1 A dx 2 A dx 3 (102) 

Therefore, the dual of d\ 4 (r) is 



*(d X4 A r ) +dX4,i(r)) = g4fcJmn(det g^y 1 ^ -detgg npV - l-— + x [(r) 



xdy* Ady j Ady k Ady 1 Ady m . (103) 

Up to first order corrections this is 



*(dXo(r) + ^Xi(r)) = *odXo(r) + 
Z' fp(r) — q{r) 
'2gZ^ 2 I Z 



^o(r)-3^i(r) +xi(r) 



* 6 d(lnr), (104) 



where *o refers to the dual in the background metric. Hence, \4 ( r ) is harmonic with respect 
to the new metric if the term in square brackets in (|104|) is zero. Thus, we set 

Xi(r) = 2^ ( Mr) ~ ^ ~ ho(r) - 3/n(r)) (105) 

Plugging in this value for Xi{ r )i one finds that the leading order corrections to the stress 
tensor coming from the cross term of dx4,i(r) and dx4,o(r), denoted by K MN , are 

K ™ = Jzk {- 4 ^( r ) J '- + w ( r ) W ™n + \p(r) ~ tfW'V)^} (106) 

for the orthogonal components and 

K v = J^W) ~ Zh^r)}^ (107) 

K °° = ~4§^ &^ ~ Zho ^ f ^ ^ 108 ) 

for the space-time components. 

The next step is to compute the linear corrections to the Einstein tensor by inserting 
these expressions for J MN and K MN into ( j5^) . The space-time components will have the 
form 

#00 - \gmR = -N (r)f(r) (109) 

Rij-^9ijR = N!(r)Vij (HO) 
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while the orthogonal components will again be of the form 



R m n-\g mn R = N 2 (r)r\r) yi !/ 



+ N 3 (r)I mn + NJr)W n 



(111) 



2" r* 

Computing the functions in ( |109|) , (|110| ) and ( |111|) is an extremely tedious exercise. We 
just state the results here. We find that 

1 



JVo(r) 



Z 3/2 



\f(r)p"(r) + (1 + 3/(r)) V(r) + (4 + 3/(r))^p(r) - |/(r) V(r) 
-(6 + 3f(r))q(r) + —f{r)h'[{r) + ±-(2 - f(r))-h[(r) - ZZ-Mr) (112) 



NAr) 



£3/2 



N 2 {r) 



ZV 2 



N 3 (r) 



z 1 / 2 



-4Z±h (r) 

\f(r)p"(r) + 2(1 + /(r)) V(r) + (6 + f(r))^p(r) 

-I(2 + /(r))V(r)-(8 + /(r))g(r) 
2 r 

+Z/(r)/i?(r) + Z(4 - 3/(r))i/ii(0 - Z(12 - 4/(r))i/i 1 (r) 
+f / W(r) + |(6 - 5/(r))~/i' (r) - Z(6 - 4/(r))i/ k) (r) 
;i + /(r)) -p'(r) + 2(2 + /(r))ip(r) - 10g(r) 

Q 7 1 1 

+ — (2 + f(r))-h[(r) - Z(6 + 3/(r))^/n(r) 

97 1 1 

+-5-/(r)-/»i(r)-3Z/(r)-j/io(r) 

2 r r 

2/(r)/f (r) + (8 + 5/(r))-/((r) + (26 + 8/(r))-U(0 



iV 4 (r) 



ZV2 



-5(1 + f{r))-q'{r) - 10(6 + /(r))g(r) 
r 

+ ^/W(r) + 10Z(2 - /(r))^(r) - 15Z(4 - f{r))^h x {r) 

Z T T 

K'V i i 

+-2"/(r)/»S(r) + 5Z(3 - 2f(r))-ti (r) - 15Z(2 - /(r))^/i (r) 
-f(r)w"(r) - -(4 + 5/(r))-iu'(r) + (6 - 2/ (r))— w(r) 



;il3) 



;ii4) 



(115) 



(116) 



2 2' ' r " r z 

Using equations (p8|)- ([l00D , (g0g )-(g0D and (IH2D-OT) we can reduce (H) to five 
coupled inhomogenous equations. We write these as 

Lo(z) = ^A(z) - \B{z) 
L 1 (z) = ^A(z)- 1 -B(z) 



12 v ' 4 
L 2 (z) = L A { Z ) + ~£(* 
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117) 



I*(z) = ~A{z) - \b{z) 

U{z) = -A{z) + B(z) 

where again z = R^/r 4 and 

n 2\rp 12 E>2 

A(z) = ^«"w 

Biz) = g2| ^[ fl2 (a(z) - ^a'(^)) 2 (1 - fa) = »W(1 - fa), (118) 

with fc = rfj/R 4 . The -l-j(z) functions are linear combinations of the metric components and 
are given by 

L (z) = 8(1 - kz)z 2 p"(z) - 2(3 - kz)zp'(z) + 3(1 - kz)p(z) 

+6(1 - kz)zq(z) - 3(3 - kz)q(z) (119) 

+24(1 - kz)z 3 h'l(z) + 12(2 - 3kz)z 2 h' 1 (z) - Qzh x (z) 
L x (z) = 8(1 - kz)z 2 p"(z) - 2(3 + kz)zp'(z) + (3 - fcz)p(z) (120) 

+2(3 - kz)zq'{z) - (9 - fcz)g(z) 

+ 16(1 - kz)z 3 h"(z) + 16(1 - 2kz)z 2 h\(z) - 4(1 + fcz)z/n(z) (121) 

+8(1 - kz)z 3 h'o{z) + 4(2 - bkz)z 2 ti {z) - 2(1 + 2kz)zh {z) 
L 2 (z) = -4(2 - kz)zp'{z) + 2(1 - kz)p(z) - 6q{z) 

-6(3 - kz)z 2 h' 1 {z) - 3(3 - kz)zh 1 (z) - 6(1 - kz)z 2 h' Q {z) - 3(1 - kz)zh (z) (122) 
L 3 (z) = 32(1 - kz)z 2 p"{z) - 32zp'{z) + 2(25 - 4A;z)p(z) 

+20(2 - kz)zq'(z) - 10(7 - kz)q(z) 

+120(1 - kz)z z b![{z) + 30(3 - TA^ 2 /^) - 15(3 + kz)zh 1 (z) (123) 

+40(1 - kz)z 3 h'o{z) + 30(1 - 3kz)z 2 h' {z) - 15(1 + kz)zh (z) 
L 4 (z) = -8(1 - fcz)zV'(z) + 8zw'(z) + 2A;zw(z). (124) 



The equations in (|117|) and (|119|) -( P-20|) look quite complicated but it turns out that they 



can be further simplified. First we note that w(z) has decoupled from the other functions 
and only appears in L4. To solve for this, write w(z) = z 1 ^ 2 w(z), where w(z) is the function 
in the inertial frame. The last equation in ( |117| ) can then be rewritten as 

™"(*) - T-V^'W - A ffi (*) = - * { U2n B \ Z \ - ( 125 ) 

I- kz Az 2 82 1 /2(1 -kz) 



Using (|118|) we see that this is the same equation as in (|79[), except the source term in (|125|) 



has an extra factor of —24. Hence the solution for w(r) is 
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w{r) = -240(r), (126) 



where 0(r) is given in (pi]). The expression in (|8l"D is quite complicated, but as we will see 



this does not matter since w(r) does not contribute to the lowest order correction to the 
entropy. 

Given ( |125|) and the arguments in the previous section, it follows that the homogeneous 
mode for w(r) corresponds to a dimension 6 operator. Consulting Table III of ]3l| 



wc sec 



that this is the scalar in the 84 representation of SO (6) with M 2 = 12. The identity of 
the two fluctuation equations (|79| ) and ( |125| ), including their sources, suggests that the two 



fields are in the same supermultiplet, although it is not clear why the sources should agree 
for T non-zero. 

Let us next consider the linear combination 

°-L Q {z) + -L t (z) + -L 2 (z) - -L 3 (z) = -A(z) - -B(z). (127) 



24 v ' 8 



This reduces to the equation 
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8(1 - kz)z 2 p"{z) - 8zp'{z) - 2(5 + kz)p(z) = — A A(z) - -B{z). (128) 

24 8 

The solution to this equation with the correct boundary conditions is surprisingly simple, 

p{z) = ~Ya a ^ z) - (129) 

If T = 0, then the homogeneous solutions for z~ l / 2 p(z) behave like z 2 and z~ l . Hence this 



mode corresponds to a dimension 8 operator. Again consulting Table III of ||31|| , this is the 
SO(Q) singlet from h a a with M 2 = 32. 

Next, we find that the first two equations and the derivative of the third equation in 



( |1 1 7|) are linearly dependent and that consistent solutions exist provided that the following 



equation is satisfied: 

(1 - kz) \2zA\z) + 6zB\z) - A{z)\ - 3(5 - 3kz)B(z) = 0. (130) 



One can easily verify that the expressions for A(z) and B(z) in ( [11 8| ) satisfy ( |130| ). 

Of course ( |130| ) must be true by general covariance; the solutions for A(z) and B(z) came 
from the equations of motion for G 5 . These were derived from the same generally covariant 
action that the Einstein equations were derived from. The fact that there are only four 
linearly independent equations means that there is an extra gauge degree of freedom. Thus, 
we can choose one of the unknown functions to be an arbitrary function. A particularly 
useful gauge choice is 

h x {z) = h {z). (131) 
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Using ( |131|) , the third equation in (|117|) reduces to 
- (2 - kz) (l2z 2 ti (z) + 6zh ) - 6q(z) = (2 - kz) (--zA'(z) + ^A(z)\ + -B(z), (132) 



where we have substituted the solution in ( |129[ ). Substituting ( |131| ) into the difference 
between the second and first equations in ( |117| ) gives 



k (lQz 3 ti (z) + 8z 2 h (z) - Az 2 q\z) + 2zq(z)) = k (^z 2 A'(z) - ^zA{z)\ . 



(133) 



Examining ( |132| ) and (|133|) , we see that the same linear combination of ho(z) and h' (z) 
appears, so we can reduce these equations to a single first order equation for q(z). Using 
( |130|) we can replace all A(z) dependent terms by B(z) and B'(z). The equation can be 
solved using standard methods, giving the solution 



, s 1 ( o/ n 2-kz r* a(z') , ; 



(134) 



where a(z) is defined in ( |118| ). An integration constant has been chosen so that q(z) has 
the correct behavior on the boundary. 

To solve for h (z) one just substitutes the solution for q(z) in ( |134| ) into ( |132j ). Using 
( |130| ) and integrating by parts, one finds 



/;„( : ) = ^- (A(z) - B{z) + 2 —^ f ^2d:J | . 



96z 



Z JO 



(135) 



Clearly, 



(PI and (pi) satisfy 



p(z)+Azh (z) -q{z) = 0. 



(136) 



Using these solutions, we can go back and find the corrections to X41 the space-time 
components of C^. From ( |95|) and (|105|) , we find that the change in X4 is 



5x4 = — / -KH ( TK A ( Z ) + \p( z ) ~ ^ z ) ~ ^ zh o( z )] ) d x ° A dx l A dx 2 A dx 3 



dz ( 1 
I r z dz' 



2gJ z 5 / 2 V12' 



24c/ J z' 5/2 



B(z') - 2 —^- f ' ^2dz") dx° A dx 1 A dx 2 A dx 3 . 



z' Jo 



(137) 



The last term can be integrated by parts. Using the relation in ( 118 ), one is left with 

1 — kz r z <r(z n 



$X4 



■dz'dx A dx 1 A dx 2 A dx 3 



(138) 



24gz 2 Jo y/z' 

In particular, we note that 8x4 = on the horizon, verifying the claim in section 2 that the 

2 




r 2 term in the probe action goes to zero at the horizon. 
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A useful check of this analysis is to explicitly compute the m 2 terms that appear in ([19]) 
at zero temperature. This term comes from corrections to C^ and the parallel components 
of the metric. Inserting these corrections into (|11]), one finds that the correction to —S/V is 



47r 2 n/i5c/ 



2 

" -.h (z) - SC 012 3(z) 



_9\fz 
Using equations ( |135|) and (|138|) , the expression in (|139|) reduces to 



(139) 



7r 2 n/i5« 



'(A {z) -B iz)+ y;°J0 dz >] 



(140) 



12g 

In the zero temperature limit we set k = 0. Using ( |118| ) and the relation 

|T 5 | 2 = 18m 2 , (141) 



one finds that the correction in ( |140|) is 



2ii 2 nu 5 a' 2 n 
m r , 

9 



(142) 



which is precisely the term required by supersymmetry. 

We close this section by listing the inhomogeneous solutions for the T = case. Setting 
k = in (H]), Q126D , ( g2§ , ( g3g ) and ([135]), we find for the gauge choice ( pl| ) 



3m 2 R 6 

w(r) = 



r 4 



gm 2 ^ 6 

p(r) = — 

m 2 R 6 . , 

lW = +^T ( 143 ) 

, . , 7m 2 R 2 

ho(r) = +-^. 



VI. ENTROPY CORRECTIONS 

In the low temperature phase the masses of the hypermultiplets are much larger than 
the temperature and therefore the contribution of the hypermultiplets to the entropy is 
exponentially suppressed. In the high temperature phase, where T » m, the masses give 
polynomial corrections to the entropy. 

At weak coupling, the entropy density S is given by 

S = ^N 2 T 3 - ^N 2 m 2 T + 0(m 4 ). (144) 
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In [29] it was shown that at strong coupling, the leading order term, So, is 



S = ^N 2 T 3 , (145) 



which is 3/4 the value of the leading term in ( |144[ ). In this section we will compute the next 
term in the series at strong coupling. 

To compute the entropy correction, one needs to compute the 0(m 2 ) correction to the 
horizon area and the temperature for a fixed Schwarzschild radius. Using the metric in 
( |101|) , we see that the 0(m 2 ) correction to the area A for a fixed r H is 

AA = IgVAg/jtAo = \"f 2 Uh(r H )^ +p(r H )\ A , (146) 

where the (x indices refer to the coordinates orthogonal to the temporal and radial directions. 
Aq is the usual AdS 5 x S$ Schwarzschild area. 

Because of the corrections to the metric, it is necessary to adjust the circumference of the 
Euclidean time circle in order to prevent a conical singularity at r = r#. The corresponding 
correction to the temperature, for fixed r#, is 



^^|(Mr,)|-?M)T. (147) 



2B?V^ y ' nj A 



Hence, using (|136|) and (|145|) , we find that the 0(m 2 ) correction to the entropy for fixed 
temperature is 

A „ [AA AT\ „ 

= 2 r f(3h (r H )^-+p(r H )^S a . (148) 



If we now use the relations in ( |118[ ), ( |129| ), ( |135| ) and ( |141[ ), we find that 



9 T(3/4) 4 m 2 ( 2 /•! x 2 [P-i/ 2 (x) - xP 1/2 (x) 
l\S - 



16 7T 3 T 2 



■■ r'Tr^ K <i49 » 



We are unaware of any analytic expression for the above integral, but the numerical result 
for AS is 

AS = -0.17UN 2 m 2 T, (150) 

which is 23% of the weak coupling result. 

This demonstrates the shrinking of the horizon area as the mass is increased. It also 
suggests that a critical temperature is reached when T ~ m. It is hoped that one can 

27 



carry out a Hawking-Page analysis of this transition |35|,|3| , by demonstrating that the free 
energy is lowered if the black-hole geometry of the bulk changes to the Polchinski-Strassler 
geometry. However, since the full 10 dimensional space is not a product space, in order to 
proceed one will require a more sophisticated approach for evaluating boundary terms in 
the Euclidean action. 

As was case for the T 3 term, the m 2 T term in the entropy is somewhat suppressed at 
strong coupling. In fact as a percentage, the suppression is more. An interesting test would 
be to carry out the analysis in |[(| for the M = 1* system. In |36|1 , the first order correction 



in l/(gN) was computed, and it was shown that the entropy increases when moving away 
from the strong coupling limit. For the M = 1* system, one should find that the absolute 
value of the m 2 T term also increases when moving away from strong coupling. 

VII. DISCUSSION 

The results presented for the backreaction required a lot of tedious calculation. But in 
the end, the results had significant simplification. This leaves us optimistic that an exact 
supergravity solution can be found. Recent progress in finding exact solutions for other 



supergravity duals |]37|j33| , |38| -^] should give one hope in succeeding in this case as well. 

It might be the case that finding an exact solution for the high temperature phase 
will be an easier task then finding the solution for the the low temperature phase. While 
the low temperature phase has unbroken M = 1 supersymmetry , it does have the extra 
complication of 5 branes wrapping 2-spheres. In any case, armed with an exact solution of 
the high temperature phase, one should be able to observe the shrinking of the black hole 
horizon down to zero area for some nonzero T = T c . Work on this problem continues. 

ACKNOWLEDGMENTS 

We thank L. Rastelli, M. Strassler and N. Warner for helpful discussions. This research 
was supported in part by the NSF under grant number PHY-97-22072 and by the U.S. 
Department of Energy under contract number DE-FC02-94ER40818. 

VIII. APPENDIX 

To get a better estimate for the temperature at which a metastable state disappears, we 
need to compute a discriminant. The derivative of (|i~9| ) is proportional to 

2Ar% 



4Ar s ( sjl - rjj/r* - 1 ) + = = + 2(r 2 - rr )(2r - r ), (151) 

rJl — rjj/r 4 



2S 



where A = (im 2 )/(2gN) » 1. Hence, to find extrema for the potential one looks for zeros 
of the polynomial 

(8A - 4)r 8 + (12 - 12A)r 7 + (-13 + 4A)r 6 + 6r 5 + (-1 - 8Ar 4 H + Ar A H )r A 

+(-12r% + 12Ar 4 H )r 3 + (13r% - AAr A H )r 2 - 6r 4 H r + r A H + A 2 r 8 H . (152) 

In fact, we want to find where two extrema coalesce, which occurs when the discriminant of 
this polynomial has a zero. The discriminant factors into a seventh order polynomial and 
and a second order polynomial in rjj, and it turns out that it is the zeros of the seventh 
order polynomial that we need. However, for large A it turns out that this polynomial is 
approximately 

2 29 A 9 (rjj) 7 - 3 3 (3351797)A 8 r 8 (r^) 5 - 2 8 3 6 ^ 6 (r^) 3 . (153) 

Hence, there is a zero at 

r H » -U 3/8 2 1/4 (3351797+ (171)(2659)v / 57) V8 A- 1/8 r , (154) 

16 v ' 



and so the solution disappears at 

Ak 9 n) 



( n 2 Y 1/8 
r H » K i ^7 r o, (155) 



where 



1 / O \ 1/8 1 la 

k = — ( — ) 3 3/8 (3351797 + (171) (2659) V^7) « 0.5522 (156) 

16 \7r 2 / v / 



In terms of the temperature and m, this is 

/ n 
\4ngNj 



( n 2 \ 3/8 
T = k\- \m\. (157) 
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